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A series of MgB2 thin films systematically disordered by neutron irradiation have been studied 
by Scanning Tunneling Spectroscopy. The c-axis orientation of the films allowed a reliable deter- 
mination of local density of state of the tt band. With increasing disorder, the conductance peak 
moves towards higher voltages and becomes lower and broader, indicating a monotonic increase of 
the TT gap and of the broadening parameter. These results are discussed in the frame of two-band 
superconductivity. 

PACS numbers: 74.70.Ad, 68.37.Ef, 61.80.Hg, 74.50.+r 



The discovery of two gap superconductivity in MgB2 [ij, |3 has renewed the interest on this mechanism that was 
theoretically predicted since the fifties. This peculiar feature occurs in MgB2 because of the presence of two sets of 
bands which cross the Fermi level. The larger gap (A^- ^ 7 meV) is related to the a bands, originated by Pxy B 
orbitals, strongly interacting with phonons, and characterized by a bidimensional Fermi surface; the smaller gap (A^r ~ 
2 meV) is related to the tt bands originated by p^ B orbitals and having an isotropic Fermi surface. Two band models 
predict that non-magnetic scattering causes pair breaking as magnetic one does in a one band superconductor; i2|| this 
is a unique effect of interband scattering, because mixing a and tt Cooper's pairs causes a complete isotropization of 
the entire Fermi surface. In the strong interband scattering limit, the two gaps should merge to one when the critical 
temperature (Tc) drops to the isotropic value of about 20-25 K.i^n^] To verify these predictions, several efforts have 
been made to evaluate the energ y g aps in samples where defects were introduced by different techniques such as 
• substitutions (Al 00,111 or C [i llE 111 El ) , irradiation [ll[l|| or in films grown naturally disordered. P5 Due to 
[ the possible simultaneous occurrence of several sources of disorder, a quantitative evidence of the role of the interband 
scattering has not yet been given. In particular, substitution may induce extrinsic effects related to charge doping, 
structural instability and inhomogeneous distribution of impurities. However, it is widely accepted that the A^^ value 
is weakly affected by disorder while A^. decreases linearly with T^,, and recently the merging of the gap has been 
observed in irradiated samples at a critical temperature (11 K) lower than the predicted one.|l4| 

In this paper, we report on Scanning Tunneling Spectroscopy (STS) measurements on c-oriented thin films where 
the defects were systematically introduced by neutron irradiation. Neutron irradiation introduces neither charge 
[ doping nor electronic or structural instabilities; furthermore, in thin films it guarantees a uniform distribution of 
defects in the entire sample. [16I] STS directly probes the Local Density of States (LDOS), whose broadening can be 
analyzed in connection with scattering mechanisms induced by irradiation. 

MgB2 thin films (2000 A thick) were grown at The Pennsylvania State University by Hybrid Physical Chemical 
Vapor Deposition (HPCVD) on 5x5 mm^ SiC substrates. Details about the samples fabrication and their basic 
characterization can be found in Ref. Neutron irradiation was carried out at the spallation neutron source 

SINQ of Paul Sherrer Institut in Villigen (Zurich). Each film was sealed under vacuum in a small quartz ampoule 
[ just after deposition; this setup allowed to perform irradiation also for long times without compromising samples 
quality. Several samples (Samples IRRIO, IRR30, IRR35, and IRR40 in the following) have been irradiated at 
[ thermal neutron fluences ranging from 6.4-10^^ to 9.5T0^* cm~^. A description of the damage mechanism, as well 
as detailed structural, transport and magnetic characterization were reported elsewhere. [T^ Low temperature STS 
experiments were performed in a cryogenic system able to operate at variable temperature and in magnetic field, 
using a Ptir home-made tip. The films were mounted on a Scanning Tunneling Microscope scanning head in inert 
helium atmosphere to preserve the surface quality. 

Table m reports the main features of the measured samples: as the fluence increases, Tc value changes from 41 K to 
17 K, and the residual resistivity po varies by two orders of magnitude. 
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TABLE L Main properties of the measured irradiated films. 



Sample 


Neutron Fluence (cm ^) 


Tc(K) 


po (fiQ-cm) 


IRRIO 


6.4-10^^ 


41.05 


1.2 


IRR30 


7.7-10^^ 


36.1 


18 


IRR35 


3.0-10^* 


22.2 


52 


1RR40 


9.5-10^* 


17.0 


87 




Figure ^ shows the typical normahzed tunneling conductance STS spectra for neutron irradiated films of different 
fluences. All the spectra were collected through a standard lock-in technique and were acquired (at T=4.2 K) by 
stabilizing the feedback loop with a tunnel current of 100 pA and a bias voltage of 20 mV. The STS spectra on 
each sample were reproducible from point to point on the sample surface, and they overlapped when normalized 
by changing the tunnel resistance. This assures that we were measuring in pure tunneling regime and guarantees a 
reduced role of possible contaminated surface layers. STS measurements were performed only when such conditions 
were achieved. A spectrum on a non irradiated film (Tc = 41 K, from Ref. ITsI) . measured in the same conditions, is 
also reported as reference. 

The plot reveals a clear monotone trend of the spectra. As Tc decreases, with increasing neutron fluences, the 
zero-bias conductance (ZBC) increases and the coherence superconductivity peaks shift to higher voltages and appear 
less pronounced. This behavior is more clearly shown in the inset to Fig.^ Only one-peak spectra were found on all 
the samples. Because in c-axis oriented films the current is injected parallel to the c-axis, the tt band contribution is 
dominant in the tunneling spectra as compared to the a band contribution. 'l^ Therefore, we identify the observed 
peaks as corresponding to the tt gap, A^. 

The increase of the ZBC is likely due to a broadening of the superconducting LDOS, as a consequence of the 
increasing disorder, which introduces subgap states and smears the LDOS divergence. The displacement of the peaks 
towards higher energies could be ascribed mainly to two effects: the LDOS broadening, supported by the simultaneous 
decrease in the peak height, and an intrinsic change of the A^r value as disorder increases and Tc decreases. 

Measurements in magnetic field perpendicular to the films surface were also performed. Figure |21 reports the ZBC 
as a function of the applied field for some samples. The plot shows that the most irradiated samples are much less 
sensitive to the magnetic field, as already observed in differently disordered samples. This result will be discussed 
in detail elsewhere. 

To get further information on the evolution of LDOS as a function of irradiation induced disorder, a quantitative 
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FIG. 2: Zero Bias Conductance as a function of the applied magnetic field (perpendicular to the surface) for samples IrrlO, 
Irr35, Irr40. 



5.0 -! 
4.5 
4.0 
3.5 
3.0 



1.5 
1.0 



0.5-1 
0.0 




□ Experimental Data 
BCS fitting curves 
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FIG. 3: STS spectra, at T = 4.2 K, on the irradiated measured MgB2 films; fitting curves, evaluated by a BCS one-gap model, 
are superimposed to the experimental data. Plots are vertically shifted for clarity. 



analysis can be developed. As a first approximation, we fitted the measured dl/dV spectra through a simple (single 
gap) BCS model with a phenomenological Dynes parameter rn.|l9j| In fact, when only the tt band gives a significant 
contribution to the tunnel current, the two-band calculation for MgB2 approaches a simple one-band BCS model. 
with and as fitting parameters. The parameter includes all the measured broadening on the superconducting 
LDOS, due to intrinsic, pair-breaking, effects, but also to extrinsic effects such as thermal noise. The results of the 
best fit curves and fitting parameters are shown in Fig. |31 

The good agreement between the fits and data in Fig. |51 confirms the adequacy of the one-band model. Although 
the high F^i value estimated for the lowest Tc sample seriously affects the quantitative reliability of the simple model, 
it furnishes anyway an indication of the film behavior. The use of a two band model, with the estimated a gap from 
the specific heat data on bulk samples irradiated through the same technique, 14] docs not change the extracted tt gap 
value substantially, but introduces more free parameters, making the overall analysis less straightforward and clear. 
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FIG. 4: A-^ values evaluated by STS spectra as a function of Tc (squares). For comparison, A,r and Ao- evaluated by specific 
heat measurements on neutron irradiated poly crystalline samples (triangles, from Ref. IT^ and circles, from Ref I3l are also 
reported. In the inset, A^/A^o vs. Tc/Tco is plotted together with a theoretical calculation. 



Figure 0] shows A^r values as a function of Tc measured from our STS spectra. A^- and Ajr evaluated by specific 
heat measurements on neutron irradiated polycrystalline samples [TM Il4j are also plotted for comparison. The A^ 
values estimated by STS are usually reported to be larger than those estimated by specific heat. This accounts for 
the slightly different A^ values obtained by the different measurements for the lightly irradiated samples. The A^ 
value estimated by STS, even considering the error bars, clearly increases as Tc decreases. In heavily irradiated thin 
films (Tc = 22 and 17 K), the extracted A^r value from the STS spectra are close to Ao- estimate d by specific heat. 

The increase of A^ is predicted by the two-band model for increasing interband scattering, [^l I2(tI| Up to now, an 
experimental proof of this effect is lacking due to the difficulty to introduce disorder in a controlled way without 
simultaneous doping effects. With the reliable gap data from the neutron irradiated samples a comparison with the 
theory is possible. In the inset to Fig.01 A^r/A^o data as a function of Tc/Tco from this work and from Refs. and 

are shown, where Atto and Tco are the tt gap and the Tc values of the unirradiated samples, respectively. The solid 
line represents the theoretical prediction of the two band model for increasing interband scattering. |2Clj For Tc/Tco 
ranging from 1 to 0.85 (ATc = TcqTc ~ 6 K), all the three data series show a similar increase of A^/A^o in agreement 
with the theoretical curve. This suggests that for low levels of disorder the main mechanism of Tc reduction is the 
pair breaking due to interband scattering. In the regime of small interband scattering, the scattering rate, Tinter: can 
be calculated as ATc/Tco = 0.2ri„(er/kBTco; with ATc = 6 K we estimate Tinter 2.5 meV. 

For Tc/Tco i 0.85 the experimental data do not agree; the A^/A^o values of this work are nearly constant and 
those from the specific heat on bulk samples decrease with decreasing Tc. In both cases, the data remain below the 
theoretical curve that continues to increase. Thus for these levels of disorder other mechanisms cause the decreasing 
of Att, and consequently the suppression of the Tc. Also for strong interband scattering a and n Cooper pairs are 
expected to mix, but the single-gap model used to extract Aj^ considers only the tt contribution. Not surprisingly, for 
the sample with Tc=17 K we found that 2A^/kBTc = 3.9, which is too high for the tt band but reasonable for the a 
band. 

In Fig. 13] the estimated Dynes broadening parameter is reported as a function of po- It increases linearly 
with Po from 0.3 to 2.8 meV. This linear correlation seems to suggest that quasiparticle relaxation processes, which 
affect Po, could enhance pair breaking mechanisms which influence Tjj. A linear increase of Tu with resistivity was 
first observed and discussed in Ref. 22;, which suggested that flnite lifetime effects due to inelastic electron-electron 
scattering vary linearly with resistivity, being enhanced by disorder. In a two gap superconductor, in principle both 
intraband and interband events can contribute to the observed scattering processes, which cannot be sorted out from 
the present data. However, the intraband relaxation rate, Tintra, can be roughly estimated by po assuming equal 
intraband relaxation rate in a and tt bands. "2^ It ranges from 0.1 to 0.6 eV (see upper scale in Fig. O, much higher 
than Yd. Interestingly F^ values are close to Tinter- For example, the sample with Tc = 36K has Tinter ~ 1.6meV 
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FIG. 5: A,r and Fd values estimated by tunnel spectra as a function or residual resistivity po (lower scale). On the upper 
x-scale, the intraband scattering Fintra has been plotted, estimated as in Ref. ^3 in the hypothesis of equal relaxation rate in 
a and n bands. The continuous line is only a guide for eyes. 

and Td = ImcV. We believe that this agreement is not fortuitous. Indeed, in a two gap superconductor the interband 
scattering is expected to cause pair breaking as scattering with magnetic impurities does in single gap superconductor; 
therefore, it is reasonable that the interband scattering can strongly influence the broadening of the tunnel LDOS of 
the TT band. Furthermore, this framework implies a linear dependence between interband and intraband scattering 
rates (suggested by the linear dependence of F^i on pQ, as discussed above) as both mechanisms scale with the defect 
density, which increases with the neutron fluence. Since the extrapolated low resistivity Tjj value is different from 
zero, part of the LDOS broadening should be ascribed to extrinsic factors such as finite noise temperature or presence 
of a contaminated surface layer. 

In conclusion, we measured, by STS, neutron irradiated MgB2 thin films grown by HPCVD. From tunnel spectra 
analysis we estimated the behavior of and F^as a function of the introduced disorder. At low level of disorder 
(Tc lowered from 40 to 36 K) A^r increases slightly as an effect of interband scattering, in quantitative agreement 
with the two band theory. Furthermore, by comparing the behavior of the Tjj parameter with resistivity we were also 
able to discuss the scattering mechanisms and reasonably infer the role of the interband scattering in determining the 
broadening of the measured density of states. 
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